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We report exact diagonalization studies of inelastic light scattering in few-electron quantum dots 
under the strong confinement regime characteristic of self-assembled dots. We apply the orthodox 
(second-order) theory for scattering due to electronic excitations, leaving for the future the con- 
sideration of higher-order effects in the formalism (phonons, for example), which seem relevant in 
the theoretical description of available experiments. Our numerical results stress the dominance of 
monopole peaks in Raman spectra and the breakdown of selection rules in open-shell dots. The 
dependence of these spectra on the number of electrons in the dot and the incident photon energy 
is explicitly shown. Qualitative comparisons are made with recent experimental results. 
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I. INTRODUCTION 

Raman spectroscopy is a powerful tool in the inves- 
tigation of the electronic properties of nanostructures-'. 
In a standard experiment, the energy of the incident 
light beam is slightly above the structure band gap (res- 
onant conditions), and the scattered light is collected in 
the backward direction (backscattering geometry). Be- 
cause of the fact that there are two photons involved in 
a Raman process, selection rules hold for the total an- 
gular momentum (or parity) of the electronic subsystem 
(A J = 0, ±2). These selection rules are different from 
the rules governing intraband absorption (A J = ±1, 
Kohn theorem), thus Raman spectroscopy allows us to 
study a different sector of the spectrum of electronic ex- 
citations in the dot. 

In the last few years, Raman mea surem ents in self- 
assembled quantum dots were 

reportecPEH. The distinc- 
tive features of these experiments are the observation of 
peaks, apparently violating the selection rule^, and the 
observation of a strong electron - LO phonon coupling 
(polaron effect^) in self-assembled dota^. The descrip- 
tion of these effects requires a higher-order theoretical 
schemd^, in which parity-violating vertices are included, 
in addition to the two electron - photon vertices. We 
notice, however, that exact calculations for Raman scat- 
tering in few-electron quantum dots, even in the lowest 
(second) order scheme, are lacking. To the best of our 
knowledge, only a few calculations for etched dots are 
available^ ^, in which the final states are properly treated, 
but the intermediate states are not. We think that a 
precise understanding of the standard Raman scattering 
processes in few-electron quantum dots is needed as a 
basis towards the description of higher-order processes. 
Thus, in the present paper we recall the orthodox 



(standard) second-order Raman scheme for light scatter- 
ing by electronic excitations in self-assembled quantum 
dots. We consider dots with up to six electrons. The 
needed wave functions to compute the Raman cross sec- 
tions are obtained by exact diagonalization in a trun- 
cated basis set of many-particle functions. Because of 
the resonant character of the process, the intermediate 
states have an additional electron-hole pair. It means 
that the largest system we should diagonalize is made up 
from seven electrons and one hole. The strong confine- 
ment regime, characteristic of self-assembled dots, makes 
it possible to reach convergence in the numerical calcu- 
lations with a relatively reduced basis set (of around 10^ 
functions). We use the Lanczos algorithm in order to 
obtain the low-energy spectrum of our Hamiltonian. 

The plan of the paper is as follows. In the next section, 
we present the model quantum dot and the way we com- 



pute the Raman transition amplitude. In section III the 
intraband excitations of the quantum dot (the final states 
in a Raman process) are described. A criterium for the 
"collective" character of a many-particle state is given, 
which is based on the multipole operators appearing in 
the off-resonance asymptotics for the Raman amplitude. 
In that section, we may appreciate how the collective and 
single-particle excitations shift in energy as the particle 
number or the confinement strength is varied. In addi- 
tion, we appreciate how, for open-shell dots, states which 
are undoubtedly charge excitations give nonzero matrix 
elements of multipole "spin" operators. Next, in section 
|IV| the interband excitations are constructed. They are 
the intermediate states in a Raman event. We compute 
the interband absorption to get an idea of the position of 
the incoming resonances in Raman scattering. In section 
|V) we present the Raman spectra. As it will be shown, 
the spectra are dominated by monopole peaks, both in 
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FIG. 1: (Color online) Schematic representation of an inelas- 
tic light scattering process. 



polarized and depolarized geometries. This is in accord 
to the fact that the dot lateral dimensions 20 nm) 
are shorter than 1/10 of the light wavelength. Partic- 
ularly interesting is the observed breakdown of the Ra- 
man selection rules in open-shell dots. Finally, in the 
last section, we qualitatively discuss the relevance of our 
calculations to the experiments detailed in Refs. |2l3l4j . 



II. THE FORMALISM 

The values of the parameters used in our model for 
self-assembled quantum dots are motivated by the ex- 
periments in Refs. |2l3l4j . For the electron mass and 
dielectric constants we took, respectively, the following 
InAs value^Sl 



Tie = 0.024 mo, and e = 14.55 cq. The in- 
plane confinement potential is assumed parabolic, with a 
characteristic frequency huje ranging in the interval be- 
tween 20 and 50 meV. This is, of course, a simplification. 
The actual confinement potential is expected to become 
fiat already at excitation energies around 100 meVi'^J^i We 
will study dots in which the number of electrons varies 
between 2 and 6. 

Under resonance conditions, the Raman scatter- 
ing trans ition am plitude is given by the following 
expressiorPSEmU^ coming from second order perturba- 
tion theory in the scattering matrix: 
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intraband excitations. Details on how the final states are 
computed, what their properties are, etc are given in the 
next section. We notice that conservation of energy leads 
to the following relation: 



Ef — Ei = hvi — hvf, 



(2) 



allowing us to express the Raman shift (r.h.s. of Eq. 
([2])) in terms of the final state excitation energy. Varying 
hvi^ a peak at a fixed Raman shift (a resonance) indi- 
cates the existence of a final state (or group of states) 
with a given excitation energy. That is how Raman spec- 
troscopy works. 

Eq. (1) contains a sum over intermediate electronic 
states. The (virtual) transitions to the intermediate 
states are caused by the electron-photon interactions. 
Because of the denominator in Eq. (jlj, when hvi is 
slightly above the dot effective band gap the sum is dom- 
inated by the resonant terms, i.e. intermediate states 
with an additional electron-hole pair, whose energies are 
Eint—Ei ~ hvi. They are described in section IV The in- 
termediate states play, of course, a role in Afi. Incoming 
and outgoing resonances, i.e. increase in Raman intensi- 
ties for precise values of hvi or hvf, are a consequence of 
resonances with given intermediate states. However, the 
position of peaks in Raman spectra is determined solely 
by the final states. Notice that our expression, Eq. ([ij, 
differs from that one used in Refs. |7I8| precisely in the 
treatment of the intermediate states, which we consider 
as many-particle and Coulomb interacting states. This 
allows us to correctly describe incoming and outgoing 
resonances in Raman scattering. 

The electron-photon vertices entering Eq. ([I]), i/^^^ 
and are single-particle operators in which matrix 

elements over photon operators were explictly computed. 
The - and -f supraindexes refer to absorption or emission 
of a photon, respectively. H'^^\ for example, is given by: 



■p|x)et/it. 



(3) 



(1) 



where we use a basis of two-dimensional oscillator states, 
cr and T, for electrons and holes qi and are the wave 
vector and polarization of the incident photon. Notice 
that the electron state in the valence band, r (conjugate 
of the hole state r), enters the matrix element in 
For H^+\ we use the relation: 



We schematically represent it in Fig. [T] hvi is the inci- 
dent laser energy. The initial state in the process is the 
ground state of the iVg-electron system, meaning that 
we are considering a process at low temperatures, where 
only the ground state is populated and hence only Stokes 
lines in the Raman spectra should be observed. The final 
states in the process, on the other hand, are quantum dot 



(4) 



where the matrix element in the r.h.s. is to be evaluated 
with the wave vector and polarization of the scattered 
photon. 



3 



J 1 i_ 

2 JNe=2 



1 -I 





J 1 L 



J 1 L 



w 12 -r' — ' — ' — ' — ' — 



CO 
CO 



4 - 





Ne=3 



r — ' — r 



10 J 



-|Ne=4 



4 J 

2 : 





T — ' — r 

J I L 



J r L 



Ik 



— I — ' — I — ' — I — ' — r 
CD 24 -H — ' — ' — ' — ' — ' — ' — ' — 

^ 18 : Ne=5 

E : 

6 - 

^ i 



J 1 L 



T — ' — r 



nn 



6 : N3=6 

4 1 
2 - 





T — ' — r 



■y 



10 2030405060 70 80 

AE^ ( meV ) 

FIG. 2: Density of energy levels in few-electron dots for 
fiue = 30 meV and excitation energy below 100 meV. See 
explanation in the text. 



III. INTRABAND EXCITATIONS IN 
SELF-ASSEMBLED DOTS 

Self-assembled InAs dots take usually the form of disks, 
of around 20 nm diameter and 5 nm high^^. We assume 
that the motion of electrons along the axis (z-direction) is 
quantized, occupying the first sub-band. For the motion 
in the perpendicular plane, we use a harmonic oscillator 
model. 

In order to compute the intraband excitations of the 
A'e-electron quantum dot we diagonalize the electronic 
Hamiltonian 
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in a basis of Slaters determinants. In Eq. E^z^ = 
h^TT^ /(2meLl), and Lz — ^ nm is the dot width. The en- 
ergy of 2D oscillator states is £„■ = 2^^. -t- |Zct-| + 1, where 
k is the radial quantum number, and I is the angular mo- 
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FIG. 3: (Color online) The low-lying intraband excitations 
in the 6-electron quantum dot. hue = 30 meV. The relative 
contribution to energy-weighted sum rules of collective states 
is given. 



mentum along the normal to the plane. The matrix ele- 
ments of Coulomb interactions, (A, /i|l/ri2|<T, x), among 
any four states of the first 20 oscillator shells were com- 
puted and stored in a file. The strength of Coulomb inter- 
actions is given by /3 = e^/(47reZe), where e is the electron 
charge, and le the electron oscillator length. With the 
exphcit values of the parameters, we get /3 = 1.75 y/htJl 
meV, where fiwe is to be written in meV also. Notice that 
the ratio between characteristic Coulomb and oscillator 
energies, P/hWf., is less than one for self-assembled dots 
(around 0.25 for huie — 50 meV). 

The dimension of the Hamiltonian matrix can be re- 
duced by restricting the basis to sectors with given values 
of the total angular momentum, L = '^^Im and total 
spin projection, Sz = Sct^^^- The dimension is fur- 
ther reduced by introducing an energy cutoff. Due to 
the strong confinement, we obtain converged eigenvalues 
of the Hamiltonian ([5| using a truncated basis of Slater 
determinants with zeroth-order (harmonic oscillator) ex- 
citation energy lower than 8 hw^. This leads to matrices 
with dimensions less than 10"* which are easily diagonal- 
ized. The algorithm used in our Fo rtran 90 code has 
many similarities with published one d^^ * "'^'' -*^. It is not 
yet parallelized, althought it may be easily adapted for 
parallel computation. 

A sample of the results is shown in Fig. |2]for dots with 
A^'e = 2 - 6, and hoj^ — 30 meV. In this figure, we present 
the density of energy levels for the excited states with 
the same L and Sz as the ground states. These are the 
final states giving the main contribution to the Raman 
cross section. Notice the ground-state degeneracies in the 
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A^e = 3 dot (quantum numbers: L = ±1, Sz = ±1/2), in 
the TVe = 4 dot (L = 0, 5^ = 0, ±1), and in the = b 
dot (L = ±1, 5^ = ±1/2). 

A more detailed view of the spectrum of excited states 
in the 6-electron dot is drawn in Fig. [3j This is a 
closed-shell system with ground-state quantum numbers 
Li = Si = 0. Using a common terminology^'^, we will re- 
fer to AL = Lf — Li = states as monopole excitations, 
AL — ±1 states as dipole excitations, AL = ±2 states as 
quadrupole excitations, etc. On the other hand, states 
with total spin variation, AS* ^ 0, will be called spin 
excitations, even if A 5*2 =0, in contrast to charge exci- 
tations which correspond to AS — 0. The cases ASz ^ 
(spin flips) will not be considered below because the Ra- 
man amplitudes for transitions to such states are very 
small^^ . Spin-flip peaks in the Raman spectra are the 
result of higher-order processes and will not be studied 
in the present paper. 

In Fig. |3] collective states are also indicated. By collec- 
tive we mean many-particle states of interacting electrons 
with signiflcant matrix elements of multipole operators. 
Let us recall the energy-weighted sum rule for charge 
excitation^i^ 



Ai?HI(/l^ilN)l' 



AL' 



(6) 

A state I/) will be conventionally called collective (a 
charge-density excitation, CDE) if Ai?/|(/|£>^^|i)p is 
greater than 5% of the r.h.s. of Eq. In contrast, a 
state with a small contribution to the sum rule is called a 
single-particle excitation (SPE). The multipole operator 
is defined as: 
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FIG. 4: (Color online) (a) Collective monopole excitations 
below 100 meV in the six-electron dot as a function of hue- 
(b) Sum rule fractions for monopole states in the few-electron 
dots and a confinement potential with hcj^ = 30 meV. The 
Kohn mode (dipole excitation) is also drawn as a reference. 



^il = E(^I^AL|M)ele^, (7) 

where the sum runs over states A and fi with the same 
spin projection^ The matrix elements {\\dAL\fJ-) are 
given elsewher^i^ESI, Making explicit the spin degrees 
of freedom, we write: 

d'aI = J2^MdAL\f^}{ey,^ + el^e^i), (8) 

in which A and /i now refer to orbital functions (no spin) . 
These operators -D^^ enter the asymptotic expression for 
the Raman amplitude in the off-resonance regim^. They 
are multiplied by the scalar product, • e*/, between the 
incident and scattered light polarization vectors. 

On the other hand, a second term in the asymptotic 
off-resonance expression^, proportional to (e^ x ef) ■ z, 
involves the operator 

D^AL = 5^(A|dAL|Ai)(el^e^T " Ai^^^^i)- (9) 



This operator may serve to distinguish collective spin ex- 
citations (spin density excitations, SDE) only when the 
gound-state total spin is Si = 0. When Si ^ 0, the corre- 
spondence is not unique, meaning that we could observe 
peaks corresponding to charge excitations in the cross- 
polarized Raman spectrum of open-shell dots. 

Notice that, in Fig. |3] the dipolar SDE's are the low- 
est states in that sector, and the monopole CDE's are 
at the top of the first group of monopole states, having 
excitation energies roughly equal to 2 hto^. 

Fig. Qa) shows the collective states of the six-electron 
dot as a function of the confinement strength, hu^. States 
have been denoted by symbols, whose sizes are propor- 
tional to their weight in the corresponding sum rule. 

On the other hand, in the lower panel of Fig. |4] the 
harmonic energy is fixed at = 30 meV, and the 
monopole states with fractions of the charge or spin sum 
rules higher than 5% are drawn as a function of the num- 
ber of electrons in the dot. There is, in addition, a re- 
markable fact in these dots, working under the strong- 
confinement regime, related to the position of the main 
CDE state: it is almost independent of the number of 
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TABLE I: The sum rule fractions represented in Fig. |4|b). 
Sl'^'"'' refer, respectively, to the charge and spin energy- 
weighted sum-rule fractions. Charge-excitation states giving 
only contributions to the spin sum rule are denoted CE's. 



electrons in the dot. Explicit values of the sum rule frac- 
tions are given in Table |T] Notice that, in the open-shell 
dots, charge excitations (collective or not) may also give 
a contribution to the "spin" sum rule. 
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a function of hcue- (b) Band gap in the few-electron dots for 
a confinement potential with huje = 30 meV. 



axis. Indeed, when — 5 nm, light hole states are 0.5 
eV higher in energy than heavy hole states. 

Under this simplification, the intermediate states are 
characterized by the total angular momentum. Lint, the 
hole spin projection, si'^\ and the total electronic spin 
projection, si'^K The main contribution to the Raman 
amplitude comes from states in which the added pair 
has angular momentum equal to zero, that is Lint — Li, 
where Li is the initial (ground) state angular momentum 
of the TVe-electron system. 

We should then diagonalize the Hamiltonian ^ with 
two additional terms 



IV. THE INTERBAND EXCITATIONS 

The intermediate states entering Eq. ([l]) for the tran- 
sition amplitude are interband excitations of the dot, i.e. 
states with an additional electron-hole pair. We will 
use a simplified description with only one (heavy) hole 

band with effective anisotropic mass, trif^ — 0.35 mg, 

^'h^^ — 0-035 mo-'^This is not a very crude hypothesis, 
which may be justified to work in the conditions of self- 
assembled dots because of the shift of light-hole states 
due to the small width of the dots along the symmetry 



fT XflfTT 

(10) 

the first one represents the single-particle energy of the 
hole, and the second one accounts for the electron-hole in- 
teractions. The basis functions are built up as products of 
Slater determinants for electrons and a single harmonic- 
oscillator function for the hole. As in the previous sec- 
tion, we define the excitation energy for these functions in 
terms of the difference of zero-order (harmonic-oscillator) 
energies. The reference is the lowest energy in the basis. 
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With a cutoff of 8 hujf, for the excitation energy (enough 
to reach convergence), matrix dimensions of around 10^ 
are obtained. The lowest 100 eigenvalues are easily com- 
puted by means of Lanczos algorithm^^B 

We show in Fig. [5] the gap renormalization effects in 
the dots as a result of varying the confinement strength 
or the number of electrons. The renormalized gap is sim- 
ply the energy difference between the lowest intermediate 
state and the ground state in the dot. Notice that the ef- 
fective gap (around 1.18 eV) is the result of adding a nom- 
inal Ega-p — 0.43 eV, the confinement energies along the 

z-direction of the added electron and hole {Ejf^ + Ez'^^), 
the in-plane confinement energies of both particles, and 
the contribution from Coulomb interactions. The depen- 
dence on hwe is almost linear, as expected, because of the 
single-particle terms in the Hamiltonian. The small red- 
shift of the effective band gap with increasing N^. , on the 
other hand, comes from the prevalence of electron-hole 
attractive interactions over phase-space filling effects in 
small dots. 

In Fig. [6] we show the interband absorption in dots 
with Nf, = 2, . . . , 6 and hujf. = 30 meV. The intention 
is to show possible incoming resonances in Raman scat- 
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six-electron dot with confinement hw^ — 30 meV. Monopo- 
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tering processes. The absorption spectrum at normal in- 
cidence is computed from the matrix elements squared, 
I (mi|iJ*^~^ |z) P, smeared out with a Lorentzian. We used 
a uniform width for the intermediate states. Tint = 0-5 
meV. Only intermediate states with the same angular 
momentum as the electronic ground states are included 
in the computation, that is Lint = Li. Two spin combi- 



nations are considered S 



and S 



~ '-'z.i -*-/ ^' '-'z.int 



q{e) I 1 /9 cW 

1/2. 



'1/2, 



Althought hijje = 30 meV is not the strongest confine- 
ment achievable in self-assembled dots, we clearly dis- 
tinguish in Fig. |6] absorption peaks arranged in shells 
separated by 2hu>e- The dispersion of peaks in the first 
shell ranges from 15 meV in the = 2 dot to a single 
(doubly degenerate) peak in the 6-electron dot. Incoming 
resonances in a Raman spectra should be then very sharp 
if experiments were conduced in a very high quality dot 
array or a single dot with the use of confocal microscopy. 

Specially interesting is the 6-electron dot, where an 
almost ideal Raman process with transitions through a 
single intermediate resonance at 1181.6 meV can be re- 
alized. 

Results for Raman cross-sections are presented in the 
next section. 



V. RAMAN SPECTRA 

The first feature of Raman spectra in self-assembled 
dots is the dominance of monopole ( AL = 0) peaks. This 
property is shared with etched dot^^^Ml^ b^t it is much 
more accentuated for self-assembled dots because of their 
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FIG. 8: Polarized Raman spectrum in the six-electron dot as 
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smaller dimensions. Indeed, the typical diameter, d, of 
InAs dots in a GaAs matrix, for example, is around 20 
nm. On the other hand, the Raman amplitude squared, 
l^/ip for a state |/) with a given AL is proportional 
to (gi sin , where (pi is the angle of incidence of 

photons, and qi ~ 0.006 nm^^ is the wavevector of pho- 
tons with energy equal to 1.18 eV. The product qid is 
roughly 0.12, meaning that only AL — peaks should 
be observed. 

We show in Fig. |7]the spectrum for a, = 6 dot and 
hiOe = 30 meV. The incident photon energy is exactly in 
resonance with the first exciton state shown in Fig. [6](e), 
that is hi/i = 1181.6 meV. The differential cross-section 
is computed from the expression: 



da ' 



E 



A 



r//(4^) 



T) + {hvf + Ef~Ei-hvif 



(11) 



where we used a Lorentzian to smear out the Dirac delta 
function expressing conservation of energy, Eq. ([2]). The 
width of final states, Fj, is assumed uniform and equal to 
0.4 meV. When the ground state is degenerate, we sum 
over all of the terms in the multiplet. 

The calculations shown in Fig. [7] correspond to the 
polarized geometry, under backscattering conditions and 
incidence angle (in vacuum) equal to 60° in order to 
strengthen multipole peaks. Contributions from AL = 
(monopole), AL = ±1 (dipole), and AL = ±2 
(quadrupole) final state excitations are included in the 
sum in Eq. (111. Dipole final states with Raman shifts 



around 30 meV (in particular, the Kohn mode) and 
around 80 meV give rise to peaks at least one order of 
magnitude smaller than the leading monopole peaks ly- 
ing in the interval from 50 to 60 meV. Quadrupole states, 
present also in this interval (see Fig. |3|, give, however, 
a negligible contribution to the cross-section. Notice 
that, under resonance conditions, peaks associated with 
single-particle monopole excitations are as strong as the 
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N =4; polarized spectra 



FIG. 9: Polarized Raman spectrum for a four-electron dot 
with confinement huje = 30 meV. The initial photon energy 
is varied in the interval where the first set of incoming reso- 
nances is expected. 



peaks corresponding to the collective states (CDE). Let 
us stress also that the main Raman peaks are shifted 50 - 
60 meV below the luminescence lines and, for this reason, 
should be easily observable. 

Polarized spectra for the six-electron dot as a function 
of the incident photon energy, hi^i, are shown in Fig. |8] 
The incident angle in this and the next figures is (pi — 
10°. Due to the special character of these processes, with 
transitions through a single intermediate state resonance, 
the dependence on hi>i is uniform, with a sharp maximum 
at the resonance. Away from the resonance, on both sides 
of it, the CDEi state gives the strongest peak. 

The Ne = 6 spectra may be contrasted with the — 4 
polarized spectra, shown in Fig. [9] The hvi intervals 
shown in both figures are the same in order to facili- 
tate comparison. Resonances with different intermediate 
states (see Fig. [6]) lead, in the present case, to a pattern 
in which the relative intensities of Raman peaks rapidly 
vary with hvi. This is the common feature of all of the 
studied dots except the = 6 one. There are also in- 
teresting facts about the Raman spin selection rules, but 
they are more evident in the next Fig. |10[ where polar- 
ized and depolarized spectra are drawn as a function of 

Monopolar spectra in Fig. [TOjare computed at incident 
photon energy exactly in resonance with Egap of each 
dot. First, we notice that for the — 2 and — 6 
(closed-shell) dots, in which the ground-state total elec- 
tronic spin is 5^ = 0, the spin selection rules deduced 
in the off-resonance regime (ORA) hold even under res- 
onance conditiond^. That is, CDE's and SPE's with 
AS* = are observed in the polarized spectrum, and 
SDE's and SPE's with AS" = 1 are observed in the depo- 
larized spectrum. This is not the case for the open-shell 
dots, in particular the = 3 and 5 dots, for which the 
polarized and depolarized spectra are very similar, and 
we can not use the selection rules in order to identify spin 



8 




hv. (meV) 



^ / f I ' 

46 48 50 52 54 56 58 60 62 
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FIG. 10: (Color online) Polarized and depolarized monopo- 
lar Raman spectra in the few-electron dots with confinement 
fttJe — 30 meV. The initial photon energy is exactly in reso- 
nance with the first intermediate state (Egap) in each dot. 



or charge excitations. In the = 4 dot, we observe the 
singlet state (a spin excitation), which can be classified 
as a SPE according to the sum rule, as a distinct low- 
energy peak in the polarized spectrum. Notice that, in 
general, the weights in the sum rule given in Table [T] are 
not indicative of the relative intensities of peaks in the 
Raman spectrum. On the other hand, the positions of 
CDE's states practically do not depend on N,., specially 
the CDE2, as mentioned in the comments to Fig. |4] and 
the SDE's peaks shift to lower energies as is increased. 

Finally, we would like to address the question about 
outgoing resonances. These resonances are defined by 
following the intensity of a fixed Raman peak as a func- 
tion of the energy of the scattered photon. In bulk sys- 
tems, outgoing resonances associated with collective final 
states were found at exactly the same energy positions of 
incoming resonances, which led to the idea that outgoing 
resonances are a consequence of a third-order process, 
in which an additional perturbation causes the decay 
of the intermediate state towards the exciton (incoming 
resonance).^'*' We would like to check what happens when 
employing exact functions in the standard (second-order) 
scheme. 
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FIG. 11: (Color online) The amplitude squared, |j4/ip, solid 
line, corresponding to the CDEi state in the = 6 dot 
and huje — 30 meV as a function of the scattered photon 
energy, huf. The absorption intensity (dashed line) is also 
shown for comparison. Notice that the x axis in the absoption 
curve corresponds to hvi. The energy of the first incoming 
resonance is indicated by an arrow. 



Fig. 11 shows the amplitude squared, , corre- 



sponding to the CDEi state, that is, the most collective 
charge-density state, in the iVe = 6 dot and Hute = 30 
meV, as a function of the scattered photon energy, hvf- 
A strong peak at hi'f « 1183 meV, 1.4 meV above the en- 
ergy of the first incoming resonance (indicated by an ar- 
row in the figure), is observed. In our scheme, this means 
a strong resonance with an intermediate state with en- 
ergy Ei^t « Ex + AE{CDEi), where Ex = 1181.6 meV 
is the position of the first absorption maximum (exci- 
ton), and AE{CDEi) = Ef{CDEi) - E,. Comparison 
with the absorption intensity, also depicted in the figure 
with a dashed line, shows that the absorption peak is 
not the strongest for this particular intermediate state. 
This means that the matrix element corresponding to vir- 
tual emission of a photon, (/|_ff(+)|mt), should be par- 
ticularly strong for this intermediate state, see Eq. ([TJ. 
We verified a similar situation with regard to the most 
collective spin-density excitation (the SDE2). The exis- 
tence of an intermediate state with such characteristics, 
energy approximately equal to the sum of two eigenener- 
gies, and a strong transition probability to the collective 
state, suggest a kind of approximate dynamical symme- 
try, and requires further research. 



VI. QUALITATIVE ANALYSIS OF REPORTED 
EXPERIMENTS 



Although there have been a number of reports of ex- 
perimental Raman studies dating from the 1990s of deep- 
etched quantum dots with large numbers (hundreds) 
of electrons (see, for example, Refs. I25l26l27p . apart 
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from tw o pre liminary theoretical studies for 6 and 12 
electrona^SES^ ^lig initial experimental investigations of 
inelastic light scattering from quantum dots containing 
few electrons commenced in 2000. In that year Chu ct 
alP observed at 5 K under resonance excitation from self- 
assembled InGaAs/GaAs quantum dots a Raman peak 
near 50 meV with a linewidth of about 25 meV. Simi- 
lar Raman spectra were observed in both polarized and 
depolarized scattering conditions without any apparent 
depolarization shift. From sample doping considerations 
and other experiments, Chu et al. deduced that there 
were about six electrons in each quantum dot and that 
the ground and first excited state of the dot were occu- 
pied. They interpreted the depolarized Raman peak as 
arising from SDEs while the polarized peak was associ- 
ated with charge density fluctuations. Our calculations 
can help clarify this assignment. As mentioned in Section 
|V]and shown in Fig. [7] monopole Raman scattering dom- 
inates in few electron dots at resonance and the Raman 
spectrum of a closed shell dot containing six electrons 
should exhibit a clear difference between SDEs in the 
depolarized scattering geometry and CDEs in polarized 
scattering. This was not the case in Chu et al.'s work. On 
the other hand, for five, but not four, electrons confined 
in the dot, the polarized and depolarized spectra are ex- 
pected to be very similar. Thus the broad 50 meV peak 
observed by Chu et al., which covers the right energy 
range according to our theory for a confinement energy 
of 30 meV, is definitely electronic in origin, is due to in- 
tersublevel transitions within the conduction band, and 
is well explained by taking five instead of six electrons 
in their dots. The lack of the expected fine structure, as 
shown in Fig. |10| in their reported spectrum is, as they 
also note, probably due to energy level variations in the 
conduction band arising from a statistical distribution of 
dot diameters in their sample. 

In 2003, Brocke et al.^^ investigated the electronic ex- 
citations in InGaAs/GaAs self assembled quantum dots 
where the dots could be filled with from 1 to 6 electrons 
by varying a gate voltage across the dots. Their resonant 
Raman scattering experiments revealed a broad Raman 
band in the energy range 40-55 meV in the polarized ge- 
ometry (there was not mention of the depolarized spec- 
trum in this paper). A feature of their results was the 
prominent peak observed at ~ 50 meV (linewidth of ^ 5 
meV) for two electrons in the dots that was assigned to a 
CDE. This peak was seen to gradually shift to lower en- 
ergy (~ 46 meV for six electrons) with increasing number 
of electrons in the dots. The shift was opposite in direc- 
tion to their expectations for CDEs and was explained 
in terms of the Coulomb interaction amongst the dot 
electrons. Apart from the Kohn mode at 50 meV, their 
energy level calculations showed the existence of other 
excitations to lower energy (in the range 40-50 meV) for 
greater than two electrons per dot, but these additional 
low energy excitations were not individually resolved in 
the experiment^. No calculated Raman spectra were 
given in this paper. Without this information, the im- 



plication of their energy level analysis performed using a 
confinement energy of 50 meV is that the clear peak they 
observed at ^ 50 meV for at least 2 electrons per dot was 
the Kohn mode. Our calculations of the Raman spectrum 
although for pure InAs dots with a confinement energy of 
30 meV produce Raman peaks in a similar energy range 
to that observed by Brocke et al., and most importantly 
show that the AL = 1 transitions including the Kohn 
mode at 30 meV give negligible contributions to the Ra- 
man intensity under resonance excitation. The present 
calculations show that in polarized scattering two CDEs 
spaced by 2-3 meV (CDEi and CDE2 shown in Fig. flOl 



are prominent, but their energies do not shift much with 
increasing number of electrons in the dot (see Fig. |4|. 
Based on these results, we can reinterpret the Raman re- 
sults of Bocke et al. as follows. Their spectra show a peak 
at ~ 45 meV and evidence that a similar peak occurs at 
^ 50 meV (this peak is more clearly seen at low and high 
electron numbers) that do not shift with electron number 
and are thus likely due to CDEs. The apparent shift of 
the ~ 50 meV peak to lower energy with increasing elec- 
tron number is an illusion created by the appearance of 
additional Raman peaks at energies between 45 and 50 
meV, especially for 4 and 5 electrons (see Fig. 10 1, that 
are unresolved in the Bocke et al. experiments. A more 
quantitative analysis than this requires further detailed 
calculations for the specific physical properties of their 
dot structure. 

More recently, resonance Raman scattering experi- 
ments have been performed on InAs/GaAs self-assembled 
quantum dots filled by n-type modulation doping with 5, 
7, and 12 electron^. Detailed results were published for 
^7 electron dots with intersublevel transitions in the ~50 
meV range. On resonance excitation a broad electronic 
Raman line was observed at 57 meV (linewidth ~15 
meV) in a near-polarized geometry that was attributed 
to intersublevel electron transitions. No detailed struc- 
ture of this Raman band was noted for ~7 electrons,*^ but 
its peak energy increased slightly with increasing num- 
ber of dot electrons from 55 mcV for ~5 electrons to 63 
meV for ^--^12 electrons, at which point some band struc- 
ture became evident .ES As opposed to the Bocke et al. 
results,!^ the apparent peak position of this Raman band 
increases slightly with the number of electrons per dot 
and thus is in accord with an assignment to CDEs. Al- 
though our calculations for InAs quantum dots extend 
only to 6 electrons they show that there are two CDEs 
expected in the vicinity of the 57 meV peak observed by 
Asian et al. and they do not shift much in energy with 
increasing number of electrons. In fact, the calculated 



polarized monopolar Raman spectrum shown in Fig. 10 
for 6 electrons exhibits a maximum at 57 meV and, when 
appropriately broadened, resembles quite well the experi- 
mentally observed Raman band for ~7 electrons. One ad- 
ditional interesting aspect of the work of Asian et al. was 
their observation of strong coupling between quantum 
dot longitudinal-optical phonons and electron intersub- 
level transitions.^ Such coupling modifies the signatures 
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of both the phonon and electronic Raman lines and thus 
for better comparison with experiment further theoreti- 
cal calculations including the effects of electron-phonon 
coupling are desirable. Once such theoretical results are 
obtained, it will become clear as to whether or not the 
broad Raman electronic lines in these samples are intrin- 
sic or, more likely, due to dot structural inhomogeneities, 
as discussed above. 

In 2005, Garcia et al.^ reported on the low energy elec- 
tronic excitations in AlGaAs/GaAs etched quantum dots 
where the dots were filled with a distribution of from 
4 to 6 electrons by modulation doping. The confine- 
ment energy was small for their sample (4 meV) and 
they observed electronic excitations in the energy range 
4-10 meV. On taking into account the difference in con- 
finement energies between their value and ours of 30 
meV (giving a 7.5x energy level scaling factor, approxi- 
mately), however, a qualitative comparison with our the- 
oretical results is still possible. Allowing for the breadth 
(linewidths of ~2-3 meV) of the two main peaks seen 
in the polarized and also the depolarized geometries, as 
also reported in subsequent work from the same group,!^ 
their results are consistent with a weighted summation 
in the two geometries of our calculated spectrum for 4, 5 
and 6 electrons giving two lower energy SDEs along with 
a nearby CDE together with a number of similar (unre- 
solved experimentally) excitations to higher energy. A 
somewhat sharper peak than the other CDEs and SDEs 
was observed at 5.5 meV in the depolarized geometry by 
Garcia et al. and it was attributed by them to an inter- 
shell monopole spin mode with AS' = — 1. Such a transi- 
tion is a special characteristic of dots containing four elec- 
trons only. Our calculations not only confirm this mode 
assignment (SDE2) but also reveal a second AS* = — 1 
mode at higher energy (SDES), as can be seen in Table 
I; these are precisely the monopole singlet (S — 0) states 
in the N = 4 dot. In our strong-confinement regime, the 
calculated SDE Raman peaks are not as prominent as, 
but have similar widths to, the main CDE Raman peaks 



(see Fig. 10 1. However, in the weak-confinement regime 



of Garcia et al., the SDE peaks could become stronger, 
because of the increasing role played by Coulomb inter- 
actions in that regime. Dot-ensemble induced inhomo- 
geneous broadening likely precluded Garcia et al. from 
observing the higher energy and weaker SDES Raman 
peak. 

In summary, our theoretical results have proved to be 
in qualitative agreement with experiment. However, in 
general, not all the detailed structure predicted in our cal- 
culations and its polarization dependence for dots filled 
with from 2 to 6 electrons has been observed experimen- 
tally, because of the wide widths (typically in the range 
5-25 meV for InGaAs/GaAs dots) of the Raman bands. 
This is partly due to present samples comprising dots 
filled with a range of electrons and also to variations in 
the size of dots in a given sample. Thus for a better com- 
parison between theory and experiment there is now a 
need for more experimental work on better defined arrays 



of dots or, better still, on single dots. Finally, we con- 
sider the obverse case of resonant Raman scattering from 
holes in self assembled dots. Such scattering has been 
only rarely investigated to date, but has been observed 
in both SiGe/Si^^ and InAs/GaAs^' quantum dots. In 
the latter work on InAs dots, p-type doping resulted in 
2-5 holes per dot depending on the as-grown dot den- 
sity. Resonant Raman scattering from such samples pro- 
duced a broad intersublevel hole excitation at ~25 meV 
(linewidth ~12 meV) for 2-S holes per dot that shifted to 
lower energies with increasing numbers of holes per dot. 
The small difference in the valence intraband energy val- 
ues obtained by Raman and photoluminescence spectra 
was explained qualitatively by the Coulomb interaction 
between electrons and holes, but the reason why the max- 
imum resonance occurs at a slightly higher energy than 
that of the hole excitation seen in Raman scattering is 
presently unknown.^ This is an area where more theo- 
retical work is needed. 



VII. CONCLUDING REMARKS 

In the present paper, we used an effective-mass descrip- 
tion of electrons and holes in self-assembled quantum 
dots, and computed, by means of exact-diagonalization 
techniques, the wave functions of initial (ground), inter- 
mediate, and final states entering the transition ampli- 
tude of an inelastic light-scattering process, Eq. ([T]). 
Polarized and depolarized Raman cross-sections in the 
backscattering geometry for dots with electron num- 
bers TVe = 2 — 6, and harmonic confinement strengths 
hio^ — 20 — 50 meV were calculated. The role of col- 
lective (charge and spin-density excitations) and single- 
particle excitations in Raman spectra was stressed. Par- 
ticularly interesting is the case of open-shell dots, where 
the spin selection rules do not hold for Raman scatter- 
ing, and the 6-electron dot, where an almost idealized 
Raman process with transition through a single inter- 
mediate resonance can be realized. We found evidence 
of approximate outgoing resonances in our second-order 
scheme, without the need of higher-order terms in the 
scattering amplitude. Existing experimental results were 
qualitatively analysed, although their proper description 
requires further work, for example, on the inclusion of 
the polaron effect in Raman scattering. 
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